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The coexistence of ferromagnetism and superconductivity in RuSr2GdCu20s was reported both 
from experiments (by Tallon et. al.) and first-principles calculations (by Pickett et. al.). Here 
we report that our first-principles full-potential linearized augmented plane wave (FLAPW) calcu- 
lations, employing the precise crystal structure with structural distortions (i.e., RuOe rotations) 
determined by neutron diffraction, demonstrate that antiferromagnetic ordering of the Ru moments 
is energetically favored over the previously proposed ferromagnetic ordering. Our results are consis- 
tent with recently performed magnetic neutron diffraction experiments (Lynn et. al). Ru t2 g states, 
which are responsible for the magnetism, have only a very small interaction with Cu e g states, which 
results in a small exchange splitting of these states. The Fermi surface, characterized by strongly 
hybridized dpa orbitals, has nesting features similar to those in the two-dimensional high T c cuprate 
superconductors. 
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I. INTRODUCTION 

The reported discovery of the coexistence of ferromagnetism and superconductivity by Tallon et. al. [[l]j|] has 
attracted a great deal of interest, since ferromagnetism breaks the degeneracy of spin-up and spin-down partners of 
Cooper pairs. || The experiments demonstrated that RuS^GdC^Os exhibits ferromagnetic ordering of Ru moments 
below a Curie temperature, Tc ~ 133 K, and becomes superconducting at a lower temperature, T c up to 40 K. From 
a theoretical point of view, Pickett et. al. Q] calculated the electronic structure of RuS^GdC^Os and discussed 
how the superconductivity, which may behave as a Fulde-Ferrell-Larkin-Ovchinniko (FFLO) type superconductor, 
can coexist with the ferromagnetism. However, the X-ray determined crystal structure used in their calculations 
was different from that recently determined by neutron diffraction [jfj - and so they did not consider the structural 
distortions (i.e., oxygen rotations) in the RuOg octahedra. 

The crystal structure of RuS^GdC^Os with a P4/mbm tetragonal space group, determined by the neutron 
experiment, || is shown in Fig. [j] (a). This is a structure similar to that of YBa2Cu307, where Y, Ba and Cu 
(chain atom) are replaced by Gd, Sr and Ru, respectively. Ru lies at a six-coordinated position in the octahedron 
composed of six neighboring oxygens (four Or u and two O apical), while Cu lies in a five-coordinated position (four 
Ocu and one O ap i ca /). Since the interatomic distance of Cu-O ap icai (2.190 A) is much larger than that of Ru-O ap i ca / 
(1.912 A), the Cu has a rather weak interaction with O ap i ca i- Thus, this may yield a two-dimensional (2D) electronic 
structure of the Cu-Ocu layer. The RuOg octahedra are rotated by about 14° around the c-axis to fill the space more 
efficiently, as shown in Fig. |l| (b). These rotations may lead to a significant reduction in the d-electron bandwidth 
due to a deviation of the angle ZRu-0_r u -Ru from 180°, and so the magnetism of Ru may be strongly perturbed 
k>( ' by the structural distortion. Hence, it is necessary to again investigate the stability of the magnetic structure in 
^ , RuSr 2 GdCu2 0s by first-principles calculations with the precise crystal structure parameters. 

With the experimentally determined structure, Q our highly precise first-principles full-potential linearized aug- 
mented plane wave(FLAPW) calculations || on three different magnetic orderings demonstrate that antiferro- 
magnetic ordering of the Ru moments is energetically favorable over the previously proposed ferromagnetic ordering. 
Furthermore, in agreement with our result, recent neutron diffraction experiments by Lynn et. al. |S| observed anti- 
ferromagnetic Ru ordering. In the present paper, we report the first- principles results on RuS^GdC^Os, and discuss 
the magnetic and electronic structures and their implications for the coexistence of superconductivity. 
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II. MODEL AND CALCULATION METHOD 



In our FLAPW calculations, we employed the experimentally determined room temperature crystal structure. @] 
For the magnetic structures, we assumed three kinds of magnetic orderings for Ru and Gd moments as depicted in 
Fig. 2: (a) a FM structure with ferromagnetic orderings for both Ru and Gd moments; (b) an AFM-I structure with 
a C-type antiferromagnetic ordering of Ru moments and a ferromagnetic ordering of Gd moments (more precisely, a 
ferrimagnetic structure); and (c) an AFM-II structure with C-type antiferromagnetic orderings for both Ru and Gd 
moments. (Note that the symmetry of both the AFM-I and AFM-II structures makes two inequivalent sublattice sites 
for Ru, O apical and Cu.) While the recent experiment || revealed antiferromagnetic alignments for both Ru and Gd 
moments with G-type orderings, in which nearest neighbors of Ru (Gd) moments along all three crystallographic axes 
are coupled antiferromagnetically, we employed the C-type antiferromagnetic ordering in order to reduce the large 
computational effort this would entail. (The G-type ordering requires doubling the unit cell of the C-type ordering.) 
This is justified because the moment alignment along the c-axis may be less important than that along the a-axis 
since the distance between the neighboring Ru (Gd) atoms in the c-axis (11.56 A) is significantly greater than that 
in the a-axis (3.84 A). Thus, the assumption of C-type magnetic structure would not alter the results or the insight 
obtained. 

Calculations were performed based on the local spin density approximation (LSDA) with the Hedin-Lundquist 
exchange correlation ||, in which the core states are treated fully relativistically and the valence states are treated 
semi-relativistically. Although the effects of electronic correlation in a strongly correlated system may be taken into 
account within a scheme such as LDA+U, RuSi^GdC^Os shows metallic character even in the RuO layer, as will be 
presented in Sec. IV, which causes that effect to be weak. Hence, while our results are restricted within LDA, they 
may be sufficient to discuss qualitative conclusions. The LAPW basis functions were used with a cut-off, \k + G\ < 3.6 
a. u., corresponding to about 1800 plane waves. Muffin-tin (MT) sphere radii in a. u. were chosen as: 2.4 (Ru), 2.5 
(Sr), 2.5 (Gd), 2.2 (Cu) and 1.2 (O). Inside the spheres, the angular momentum expansion was truncated at I = 8 
for wave functions, charge density and potential. Integrations were performed over a first Brilliouin zone using 126 
special fc-points, corresponding to 16, 16 and 32 fc-points in the irreducible Brilliouin zone for the FM, AFM-I and 
AFM-II structures, respectively. 



III. RESULTS 



A. Magnetic Moments 

The calculated magnetic moment of each atom inside the MT sphere for the FM, AFM-I and AFM-II structures is 
given in Table Q. The moment of Gd is close to 7 [is - as expected. In the AFM-I structure, the magnitude of the Ru, 
Oapicai and Cu magnetic moments with opposite spin directions is slightly different since the ferromagnetically ordered 
Gd moments break the symmetry of the two sublattice sites. Among all three structures, the Ru magnetic moments 
are almost of the same magnitude, about 1.55 hb- Note that since the magnitude depends on the MT sphere radius, 
the rather large MT sphere radius (2.4 a.u.) we used gives rather large values. However, the experimental value 
(1.18 iib) of the Ru magnetic moment determined by neutron diffraction [Bj is smaller than our results. This may be 
attributed to possible canted antiferromagnetic ordering of the Ru magnetic moments, for which further investigations 
are necessary. The induced magnetic moments of O ap icai and Or u are very sensitive to the Ru magnetic structures, 
while the Cu and Ocu moments are negligibly small. 



B. Total energy 

Total energies for all three FM, AFM-I and AFM-II structures are also given in Table |. From a comparison 
between the FM and AFM-I structures, we found that an antiferomagnetic ordering of Ru moments is favored over the 
ferromagnetic ordering. The total energy difference, 22.5 meV/Ru-atom, between both FM and AFM-I corresponds 
to the Neel temperature (Tjv) of 73 K, crudely estimated by mean field theory with S = 3/2. Although the calculated 
Tjy is lower than the experimentally determined TV, 136 K, ]8| our prediction for the magnetic structure is consistent 
with the experiment by Lynn et. al. Q Further, we found that an antiferromagnetic ordering of the Gd moments 
(AFM-II) reduces the total energy, but only by 2.3 meV/Gd-atom, from that of the AFM-I structure, corresponding 
Tjv « 7 K (assumed S — 7/2), which is approximately in agreement with the observed low Neel temperature, 2.5 K, 
H for the Gd moment alignment. 
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IV. DISCUSSION 



A. Density of states 

Figure § shows the total density of states (DOS) and partial DOS of each atom for the FM, AFM-I and AFM-II 
structures. The peaks at 2 eV above the Fermi level (Ep) and 3 eV below Ep in the total DOS for all three structures 
correspond to Gd 4/ states which are strongly localized, and their states do not significantly affect the electronic 
structure around Ep. An exchange splitting in the total DOS in the AFM-I structure is due to the ferromagnetism 
of Gd, but the difference between majority and minority spin DOS is very small except in the region of the localized 
Gd 4/ states. Note that the total DOS for all three structures shows metallic character. 

From the partial DOS in Fig. ^, we can see that the magnetism of Ru is dominated by antibonding t2 g states, and 
the major Or u and O ap icai peaks correlate strongly with those of Ru. In contrast, the partial DOS for Cu d and Ocu 
p states, which are strongly hybridized, are insensitive to the magnetism of Ru. The exchange splitting of the Ru d 
states for all three structures is about 1 eV, while the induced exchange splitting in the Cu d states is smaller by two 
orders of magnitude - 0.01 to zero eV. Our results for the exchange splitting and magnetic moments of the Cu d states 
is smaller, even in the FM case, than those in the undistorted FM case reported by Weht et. al. |l0| This may be due 
to the differernt structure parameters employed in the two calculations. The fact that the Ru-Cu distance, 4.10 A, in 
our calculations is larger than in their case, 3.58 A, results in a rather weak magnetic interaction and leads to smaller 
exchange splittings and magnetic moments in the Cu d states. Futher, the relative atomic position of O ap i ca i between 
the Ru and Cu may be sensitive to the magnetism. |lC| ] However, the exchange splitting and the magnetic moments in 
the Cu d are quite small. The small exchange coupling is due to the unique electronic structure of the layered Ru t2 g 
and Cu e g states separated by p orbitals of O ap i ca i, as discussed by Pickett et. al., j| and is roughly valid regardless 
of the magnetic ordering of Ru moments: Ru ti g states couple only to the p x and p y orbitals of O ap i ca i, but do not 
couple to the Cu-0 d x 2_ y 2 states. 

B. Magnetism of Ru 

We focus here on the magnetism of Ru. Figure ^ shows the projected DOS in real space for Ru d, Or u p and O ap i ca i 
p states for the FM, AFM-I and AFM-II structures. The x and z-axes for the Ru (Or u ) DOS are chosen as directions 
to the neighboring Or u (Ru) site and c-axis, respectively. Compared to the DOS of the FM structure presented by 
Weht et. al., JlQl which did not include a structural distortion due to RuOe rotations, our calculated majority d xy 
states (of., Fig7~|] (a)) are more localized below Ep, due to deviation of the angle ZRu-Oij u -Ru from 180° and the 
elongation of Ru-Ofl u bonds by the RuOe rotations. This results in the d xy states being almost fully occupied, and 
in a depletion of the majority spin hybridization channel at Ep. In the AFM-I and AFM-II structures, the majority 
d xy and 4rz(yz) states on a Ru site can hybridize with the minority spin states on neighboring Ru sites through On u 
p y and Or u p z orbitals, respectively, by the superexchange mechanism. The d xz ( yz ) states also couple strongly to 
the Qa P ical Px(y) orbitals. The charge configuration of Ru seems to be t\ g with a high spin state, namely close to 
Ru 5+ . Generally, the superexchange interaction between fully occupied t\ g states in nearest neighbor magnetic ions 
tends to be stabilized with an antiferromagneticic ordering of their moments. On the other hand, the double exchange 
mechanism through itinerant electrons may not be favorable, since the majority t2 g states are strongly localized and 
almost fully occupied - as seen in Fig. || (b) and (c). Therefore, the magnetism of Ru in RuSr 2 GdCu 2 08 is determined 
by the superexchange mechanism. 

C. Band structure and Fermi surface 

The calculated band structures for RuSi^GdC^Os are shown along high symmetry directions in Fig. [s]. The 
majority and minority spin states in AFM-II are degenerate due to the symmetry of the antiferromagnetic Ru and Gd 
moment alignments, as shown in Fig. ^. In the AFM-I structure, a small difference between majority and minority 
spin states was observed due to the assumed ferromagnetic Gd moments, but its difference is negligibly small. The 
band structure of AFM-I is practically the same as that observed in the AFM-II structure. The bands around Ep for 
all three structures arise from the Ru-O^-Oapica; bands (solid light circles) and Cu-Ocu bands (solid dark circles). 
In the AFM-I and AFM-II structures, there are three bands crossing Ep - as indicated in Fig. |^ (c). They arise 
from a Ru-Oij„ band with antibonding d xy — p y (dp?:) orbitals, and two Cu-Oc« bands composed of antibonding 
d x 2_ y 2 —p x (dpa) obitals with even and odd symmetry derived by a mirror operation in the Cu-Oc« bilayer. The two 
Cu-Octi bands are almost the same as those in the FM structure. As expected from the structural similarity between 
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RuSr 2 GdCu208 and YBa2Cu30y, their band structures arising from Cu-Ocu layers are close to each other, showing 
a strong two-dimensionality. The energy states are little changed upon going from T to Z along the z direction. 
Also, the strong dpa hybridization leads to a wide bandwidth of 10.1, 9.7 and 9.7 eV for FM, AFM-I and AFM-II, 
respectively. These values are almost the same as observed in the 2D dpa bands in YBa 2 Cu 3 07, 9.0 eV. 0|| The 
slightly larger bandwidth may be due to the smaller interatomic distance (1.93 A) of Cu-Ocu bonds and the smaller 
buckling angle (5.7°) in the CuC>2 plane of RuSi^GdC^Os, compared with those of YBa2Cu307 (1.95 A and 7.8°). 

From a comparison of band structures between the FM and AFM-I structures, we found that the antibonding dpa 
bands are not influenced significantly by the magnetism of Ru except around the M point in the zone. Around M, 
where they possess small components of Cu d 3z 2_ r 2 character, the Cu e g states couple to the Ru e g states through 
the O a pi ca i p z orbitals, but do not couple to Ru t 2g states. Therefore, the quite small exchange splitting in the dpa 
states may be mediated through a more indirect coupling path with O ap i ca i-0 r u . Note that the Cu t 2g states are 
coupled to the Ru t 2g states through the O ap i ca i p x (y) orbitals, but the Cu t 2g states are almost fully occupied, so the 
hybridization channel at Ep is depleted. 

The calculated Fermi surfaces (FS) at k z = for the AFM-II structure are shown in Fig. || (a). For convenience, 
a sketch of the FS is given in Fig. | (b), which makes for an easier comparison with the high T c results. Two closed 
FS centered at the zone corner arising from the Cu-Ocu bands and a closed FS from the Ru-0/j u band, are found. 
The inner and outer FS (1 and 2 in Fig. ^| (b)) correspond to antibonding dpa states with even and odd symmetry 
under z reflection, respectively. Note that the FS has no exchange splitting, due to the antiferromagnetic ordering of 
the Ru and Gd moments, implies that the superconductivity need no longer be of the FFLO type - in contrast to the 
FM case. Q A folding in fc-space due to the antiferromagnetic ordering of Ru makes intersections between the bands 
(a, b and c in Fig. ^ (b)), showing the band-band interaction. As observed in the undistorted FM, B the FS clearly 
demonstrates similar nesting structures with those in the other high T c cuprate superconductors |0,0]; these will 
give rise to singularities in the generalized susceptibility, and lead to possible anomalous behavior of the electronic 
properties. Further, the FS (3 in Fig. ^ (b)), derived from Ru d xy states, shows a pocket of electrons centered at T, 
which may lead to the superconductivity in the Cu-Ocu layer with hole-doped character, as observed in experiments 

Si 

In contrast, the FS in the FM structure is complicated due to the Ru-O/jn bands. The majority spin Ru d xz r yz } 
bands (4 in Fig. || (a)) become more dispersive and cross at Ep compared to the AFM-I structure, which makes an 
additional hole pocket centered at T, leading to be a less hole-doped character in Cu-Ocu layers. Hence, the FM 
structure may show low T c or non-superconductivity. 

V. SUMMARY 

We have performed first-principles FLAPW calculations to investigate magnetic structures in superconducting 
RuSr2GdCu2 08- Contrary to the previously ferromagnetic RuSi^GdC^Os structure proposed from susceptibility 
experiments jl],|2| and first-principles calculations, 0] we found that antiferromagnetic ordering of the Ru moments is 
energetically favored over ferromagnetic ordering. Our results are consistent with those of a recent neutron diffraction 
experiment, ^| in which the antiferromagnetic ordering of Ru moments was predominantly observed. The magnetism 
arises from Ru t2 g states through the superexchange mechanism. The Ru t 2g states have a very small interaction 
with the Cu e g states - which results in quite a small exchange splitting in the antibonding dpa states. The Fermi 
surfaces derived from the antibonding dpa states have similar nesting structures with those in the 2D high T c cuprate 
superconductors, which may give a basis for superconductivity in this material. Further experimental and theoretical 
investigations are necessary to confirm this conclusion. In any case, the coexistence of antiferromagnetism and 
superconductivity in RuS^GdC^Os is clearly allowed. 
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TABLE I. Total energy difference, AE (in meV/cell), and calculated magnetic moments, m (in /is), in the MT sphere of 
each atom for the FM, AFM-I and AFM-II structures of RuSr 2 GdCu208- The symmetry of the AFM-I and AFM-II structures 
makes two inequivalent sublattice sites for Ru, O ap icai and Cu. 
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FIG. 1. (a) Crystal structure of RuSr2GdCu20s with P4/mbm space group; (b) Top view of Ru-Or u layer. 

FIG. 2. Magnetic ordering of Ru and Gd magnetic moments for (a) FM, (b) AFM-I and (c) AFM-II structures of 
RuSr2GdCu208, where the Cu, Or u , Ocu and O ap icai atoms are not given. 

FIG. 3. Total density of states (DOS) and partial DOS of Ru d, Or u p, Oapicai p, Cu d and Ocu p orbitals for the 
FM, AFM-I and AFM-II structures of RuSr2GdCu208- Solid and dashed lines represent majority and minority spin states, 
respectively. A vertical dotted line denotes the Fermi level. 

FIG. 4. Projected density of states (DOS) in real space for Ru d, Or u p and O ap icai P orbitals for the FM, AFM-I and 
AFM-II structures of RuSr2GdCu20s. The x -axis for the Ru d (Or u p) DOS are chosen as directions to neighboring Or u 
(Ru) sites. 

FIG. 5. Calculated band structure along high symmetry directions for the FM, AFM-I and AFM-II structures of 
RuSr2GdCu208- The bands originating mainly in the Cu-Ocu layers are represented with solid dark circles. 

FIG. 6. (a) Calculated Fermi surfaces (FS) of RuSr2GdCu20s for the AFM-II structure at k z = and (b) its schematic 
illustration, which makes for an easier comparison with the high T c results. 
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Fig. 2 K. Nakamuira, et. al. 
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Fig. 3 K. Nakamura et. al. 




I I I I I I I I I I I I I I I 

-2-1 1 2-2-1 1 2-2-1 1 2 



Energy (eV) Energy (eV) Energy (eV) 

(a) FM (b) AFM-I (c) AFM-II 



Fig. 4 K. Nakamura, et. al. 
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Fig. 5 K. Nakamura, et. al. 
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